Astronomy & Astrophysics manuscript no. Fg071 


February 2, 2008 


(DOI: will be inserted by hand later) 





Line shape diagnostics of Galactic 26 AI 

K. Kretschmer 1 , R. Diehl 1 , and D. H. Hartmann 2 

1 Max-Planck-Institut fur extraterrestrische Physik, Postfach 1312, D-85741 Garching, Germany 

2 Department of Physics & Astronomy, Clemson University, Clemson, SC 29634-0978, USA 

Received 4 July 2003 / Accepted 6 November 2003 

Abstract. The shape of the gamma-ray line from radioactive 26 A1, at 1808.7 keV energy in the frame of the decaying isotope, 
is determined by its kinematics when it decays, typically 10 6 y after its ejection into the interstellar medium from its nucleosyn- 
thesis source. Three measurements of the line width exist: HEAO-C's 1982 value of (0 + 3) keV FWHM, the GRIS 1996 value 
of (5.4 ± 1.3) keV FWHM, and the recent RHESSI value of (2.0 ± 0.8) keV FWHM, suggesting either "cold", "hot", or "warm" 
26 A1 in the ISM. We model the line width as expected from Galactic rotation, expanding supernova ejecta, and/or Wolf-Rayet 
winds, and predict a value below 1 keV (FWHM) with plausible assumptions about 26 A1 initial velocities and expansion history. 
Even though the recent RHESSI measurement reduces the need to explain a broad line corresponding to 540 km s -1 mean 
2< A1 velocity through extreme assumptions about grain transport of 2< A1 or huge interstellar cavities, our results suggest that 
standard 26 A1 ejection models produce a line on the narrow side of what is observed by RHESSI and INTEGRAL. Improved 
INTEGRAL and RHESSI spatially-resolved line width measurements should help to disentangle the effects of Galactic rotation 
from the ISM trajectories of 26 A1. 

Key words, nuclear reactions, nucleosynthesis, abundances - gamma rays: observations - supernovae: general - ISM: super- 
nova remnants - stars: formation. 



1. Introduction 

Current Galactic nucleosynthesis reveals itself through the de- 
cay of 26 A1, one of its radioactive by-products with a mean life- 
time of 1.04 x 10 6 y. 26 A1 undergoes /3 + -decay into an excited 
state of 26 Mg, which de-excites through emission of a gamma- 
ray photon at 1808.7 keV. This gamma-ray li ne has been ob- 
served and image d throu ghout the Galax y IM ahonev e i tall 
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1982 | iDiehletail Il995t loberlackl Il997t Iknodlsed eretaTI 
1999; Pluschk e et al.Ll200ll) . Sources of 26 A1 may be AGB stars 



and novae, but massive stars (via core-collapse supernovae 
and winds from Wolf-Rayet stars) have been foun d the most 
plaus ible and probably dominating sources JPrantzos & DiehlL 
Il996l) . The rather irregular 26 A1 emission along the plane of 
the Galaxy, and its consistency with the patterns of tracers of 
massive-star activity, is t he main argument for favouring mas- 
sive s t ars as the sourc es (Diehl etal., 1996; Knodl seder et all 
199 ( 



Ipltischke et al. . 2001]). Flux measurements have been 

employed to study the nature of the sources, comparing with 
predicted 26 A1 yields from models of the source types. The 
amount of 26 A1 present in the ISM of the Galaxy has been es- 
timated at ss 2 M , and used to argue for the roles of different 
source types. But the uncertainties about the spatial distribution 
and total number of nucleosynthesis events add to source yield 
uncertainties, providi ng only qualitative argum ents for the na- 
ture of 26 A1 sources (Prantzos & Diehl, 1996). Therefore, lo- 



Send offprint requests to: K. Kretschmer, e-mail: kkr@mpe.mpg . de 



cally constrained candidate source populati ons have been stud- 
ied, such as in the C ygnus region (Knodlsed er et all [2000; 
Pluschke et all 120011) . In such a case, the distance to the 
sources is constrained to a smaller interval, along with the ra- 
dial velocity due to Galactic rotation. 

Measurements of Galactic 26 A1 with high-resolution spec- 
trometers have produced somewhat controversial results: The 
initial discovery with the space-borne HEAO-C Ge spectrom- 
eter had rep orted a narrow line ( intrinsic width (FWHM) 
(0 + 3) keV; Ma honev et all Il984l) . But the GRIS balloon- 
borne Ge spectrometer found the line to be significantl y broa d- 
ened (intrinsic FWHM (5.4 ± 1.4) keV: iNava etaP Il996)). 
Such a broad lin e, however, cannot be explained easily (see 
Che n et al l ll997l) . If the origin of the line broadening was ther- 
mal, the 26 A1 decay region would need to be at a temperature of 
=s 4.5 x 10 K. Alternatively, 26 A1 isotopes would have to main- 
tain a mean velocity around 540 km s _1 over their 1 My decay 
time, travelling kpc distances at those speeds. One may either 
assume that a substantial fraction of 26 A1 is injected into such 
rather large interstellar cavities, or that a substantial fraction of 
26 A1 condenses onto dust grains before deceleration, so it main- 
tains its momentum throughout passages of supernova remnant 
shells or other obstacles. The velocity could also be a result 
of re-acceleration of dust grains by interstellar shocks in the 
neighbourhood of the source, allowing them to maintain a high 
velocity over the 26 A1 decay time scale ( Stu rner & Naval ll999). 
None of these explanations is straightforward or without prob- 
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lems. A firm measurement of the Al gamma-ray line width 
is desirable, before questioning our understanding of 26 A1 fate 
from its production sites until decay in interstellar space. 

New measurements have been obtained recently with Ge 
spectrometers aboard the hi gh energy solar spectroscopy im- 
ager RHESSI dSmifhll2003l) and the INTEGRAL observatory. 
The RHESSI result was derived through Earth occultation anal- 
ysis of data while pointing at the sun. They obtain an interme- 
diate intrinsic line width of about (2.0 + 0.8) keV FWHM. The 
preliminary INTEGRAL measurement also suggest s a narrow 
line, but systematic uncertainties are still large dPiehl et all 
120031) . Within the given uncertainties, the current set of mea- 
surements is mildly inconsistent. But more measurements have 
been recorded with INTEGRAL'S SPI Ge spectrometer already, 
so that longitude-resolved line width results are at the horizon. 

We explore the potential of 26 A1 decay line shape mea- 
surements for diagnostics on the sources of 26 A1 and their in- 
terstellar environment. We follow up on earlier analysis by 
iGehrels & Chenl lll996l) . who showed that the structure of the 
Galaxy should be reflected in the position of the line due to 
Doppler shift from Galactic rotation. In this paper we em- 
ploy an updated spatial distribution model for the plausible 
sources of 26 A1 in the Galaxy. We also account for the fate of 
ejected 26 A1 in interstellar space through different model vari- 
ants which reflect ejections by winds and/or supernovae into 
cavities of plausible sizes for the massive-star environment of 
the 26 A1 sources. Our aim is to illustrate the diagnostic power of 
measuring the 26 A1 gamma-ray line's position and width; this 
should be feasible with current instrumentation through imag- 
ing spectroscopy, even if a detailed decomposition of the line 
shape may still be beyond reach. 

2. 26 A1 Sources in the Galaxy 

Our model for the distribution of 26 A1 in the Galaxy is based 
on two separate aspects: the distribution of nucleosynthesis 
events in the Galaxy, and the distribution of 26 A1 in interstellar 
space following an individual nucleosynthesis event. For each 
of these, spatial as well as velocity-space densities have to be 
considered. 

The angular distribution of 26 A1 1809 keV emission on 
the sky correlates well with tracers of ionisation, such as H a 
or free-free emission. Therefore we adopt a three-dimensional 
model for the space density of free electrons as our parent dis- 
tribution for 26 A 1 sources in the Galaxy . Such a model has 
been derived by iTavlor & Cordesl 0993I) from pulsar disper- 
sion measure observations (a cut through this model along 
the Galactic plane can be seen in Fig. ??); it was updated in 
ICordes & Laziol (l2002). Alter native spatial models , such as the 
smooth, axisymmetric one bv lGomez et alJ J200lh will proba- 
bly yield similar line shape results. 

The Doppler shifts due to Galactic rotation can then be de- 
termined from the Galactic rota tion curve. For our mod el we 
used the results obtained by lOlling & Merrifieldl ll200Cl) from 
fitting radial velocity measurements with a five-component 
mass model of the Galaxy, consisting of a stellar bulge, a stel- 
lar disc, two gas discs (H i, H2) and a dark-matter halo. Because 
their determination of the rotation curve depends on the values 
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Fig. 1. M odel of the free electron density in the galactic plane 
llTavlor & Cordesl 1 19931) . used as 26 A1 source density distribu- 
tion. Dotted lines represent -60°, -30°, -4°, 0°, 4°, 30°and 60° 
galactic longitude. Hatched areas illustrate the eastern/western 
parts of the inner galactic region. 

of the distance to the galactic centre Rq and the local circular 
speed ©0, which are not known to a high degree of precision, 
they allowed these parameters to vary over a broad range of 
values. In view of these uncertainties and because we are inter- 
ested mainly in the inner Galaxy, we approximated the rotation 
curve given by Oiling & Merrifield for the IAU standard values 
= 8.5 kpc and ©0 = 220 km s with the radial dependence: 

\v\(R) = 220 km s- 1 ■ [1 - exp(-fl/l kpc)] 

the velocity vector being parallel to the plane and perpendicular 
to the vector pointing from the galactic centre to the source 
location. 

Superimposed on Galactic rotation is the motion of freshly 
synthesised radioactive material due to the parental supernova 
explosion or the ejecting Wolf-Rayet star wind and its slowed- 
down motion in the ISM before decay, i.e. within w 10 6 yr. 
Concentrating first on supernovae, we adopt a particular expan- 
sion behaviou r: Recent hydrodynam ic simulations of type II 
supernovae by Kifonidis et al. ( 2003) find that the expansion of 
the bulk of nucleosynthetic SN products such as 26 A1 may be 
at velocities less than 1200 km s . To reflect their results, we 
allow 26 A1 to expand freely with a velocity of 1500 km s until 
it reaches the radius of the SN reverse shock formed by circum- 
stellar interaction. After this point, we expand 26 A1 at the veloc- 
ity of the blast wave shock; this gives us a conservative estimate 
because 26 A1 is likely to move slower than the forward shock. 
For our model of SNR dynamics from circumstellar interaction, 
we adopt the values for Kepler's sup ernova remnant shoc k po- 
sitions and velocities given bv lMcKee & TrueloveNl995l) . 

At present, our model does not include other 26 A1 sources; 
Type Ib/Ic supernovae and Wolf-Rayet stars ( Pr antzos & D iehl , 
1996), which eject matter at similar or even higher speeds than 
Type II SNe dMellema & LundavisdEoollPrinia et ailll990t 
Gar cia-Segura et all Il996l) . Also, the interaction of ejected 
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matter with the surrounding medium depends on the star for- 
mation history of the source region, where bubbles forming 
around groups of young massive stars play a potentially large 
role. Cavities exte nding over several hu ndred pc have been ob- 
ser ved in galaxies I Oev & Clarkel! 19961) . and the Eridanus cav- sg L oi 
ity feurrows et all 19931) presents us with a nearby example of 
such a cavity, extending from the Orion star forming region 
to very near the Sun. Matter ejected into such a low-density 
bubble could expand almost freely until reaching the bound- 
ary whereas typically assumed ISM densities («h - 1 crrT 3 ) 
would slow it down quite rapidly. To keep the model simple at 
first, we adopt above SNR model as typical for 26 A1 sources; 
refinements will be discussed in our subsequent studies. 

With these assumptions, the ejected 26 A1 moves freely for 
w 2 kyr, then decelerates with the SNR's shell. The shell 
reaches a radius where it dissolves in the ISM at an age com- 
parable to the lifetime of 26 A1 , when a significant fraction has 
therefore already decayed. We note that the expansion velocity 
drops below the characteristic rotational velocity of 220 km s _1 
at ss 40 kyr, when 96% of the 26 A1 is still left. Therefore the 
contribution from expansion to the overall line width will be 
rather small in our model. 
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Fig. 2. 



26 A1 line intensity as a function of Galactic longitude 
and gamma-ray photon energy, with the source longitude pro- 
file superimposed 
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3. Line Shape Diagnostics 

We obtain simulated sky maps and spectra from a Monte Carlo 
scheme: 26 A1 source locations are chosen randomly from a spa- 
tial distribution proportional to the free electron density, dis- 
tributing candidate source positions within a volume centered 
on the Galaxy and extending 25 kpc in the plane and 7.5 kpc 
perpendicular to the p lane. The free-electron d ensity within this 
volume is taken from T avlor & CordesHl993l) . For each nucle- 
osynthesis event, a random age is chosen within the interval 
[0, 10 7 yr]. From this we evaluate the intrinsic velocity distri- 
bution and size of the 26 A1 source, following the above expan- 
sion model. The age of the nucleosynthesis event thus deter- 
mines extent and intrinsic velocity of its ejecta, as well as their 
1.809 MeV luminosity. We represent each event by 2 10 mass 
elements to reflect its spatial extent. 

With the observer at Rq = 8.5 kpc and 0n = 220 km s _1 , 
we obtain viewing direction and radial velocity of 26 A1 sources. 
Direction and radial velocity give us the coordinates of the 26 A1 
source mass element in a data space of 26 A1 decay luminosity 
as a function of longitude, latitude and photon energy. 

Projecting this data volume onto the longitude-energy 
plane, we obtain a map of 26 A1 line intensities as a function of 
Galactic longitude and observed photon energy (See Fig. ??). 
The intrinsic velocity spreads and spatial source distributions 
lead to a line blurring of sa 0.3 - 0.5 keV along the plane of the 
Galaxy. Galactic rotation leads to the symmetric line shifts east 
and west of the Galactic centre. Deviating from the trend of 
decreasing surface brightness with increasing separation from 
the Galactic centre, we notice a prominent flux enhancement 
around -70°, which corresponds to the direction tangential to 
the Carina-Sagittarius spiral arm, which is also clearly visible 
in the longitude profile of our source distribution model, super- 
imposed as a histogram in the upper half of the figure; see also 
Fig.??. 
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Fig. 3. Spectrum of the inner galaxy (/ € [-30°, 30°]. The 
full line is the source spectrum, the dotted line the best-fitting 
Gaussian, and the dashed line results from a convolution with a 
2.8 keV FWHM Gaussian adopted for instrumental resolution 
of the measuring detector, b e [-5°, 5°]) 



In order to compare with a measured line profile, we in- 
tegrate our simulated skymap of line energies and intensities 
over a region of interest of and obtain a resulting spectrum 
corresponding to an observation of this region without spa- 
tial, but with perfect spectral resolution. The spectrum shown 
in Fig. ?? represents a region-of-intere st chosen to re flect the 
RHESSI analysis of the inner Galaxy JSmitrl l2003h . Fitting 
a Gaussian to this spectrum yields an equivalent line width of 
1.0 keV (FWHM). The deviation from a Gaussian shape, which 
appears clear in our model spectrum (Fig. ??, continuous line 
versus dotted line) is probably too small to be detectable with a 
realistic Ge detector of instrumental resolution 2.8 keV FWHM 
(approximate for SPI on INTEGRAL, dete rmined by in-orbit 
measurements of detector background lines; Atti e et all 120031) . 
convolution with the instrument response suppresses the differ- 
ence between the measured flux and a Gaussian by a factor of 
~ 2000 (see dashed line in Fig. ??). 
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Fig. 4. Illustrations of lines from different Galactic regions: The 
histograms (dotted/dashed lines) show the spectra of the east- 
ern/western part of the inner galactic region (i.e. ±5° in lati- 
tude, +[4° - 30°] in longitude) Energy -resolution limited spec- 
tra (obtained by convolving with 2.8 keV FWHM Gaussians) 
are shown as dash-dotted/long-dashed lines. 



Nevertheless, spatially resolved observations can - at least 
partially - separate the rotational effect from the total line 
broadening. Fig. ?? compares the model spectra of the east and 
west part of the inner galactic region, (defined by b e [-5°, 5°] 
and I 6 [4°, 30°], / € [-30°, -4°], respectively) If we assume 
an instrument with a spectral resolution of 2.8 keV FWHM, 
the observed spectra are again almost indistinguishable from 
Gaussians, but we expect their centroids to be shifted by a sig- 
nificant energy offset. For our choice of parameters, she shifts 
amount to -0.25 keV for the eastern and +0.25 keV for the 
western part. A real measurement would yield a line width that 
is lower than the width obtained for the entire inner Galaxy. 
When we fit Gaussians to these resolution-limited spectra, we 
obtain as 0.9 ke V, compared to « 1 .0 keV for the whole inner 
region (b e [-5°,5°], * € [-30°, 30°], Fig. ??). This is due 
to the fact that partial-region spectra only include a lower dy- 
namic range of velocities. 

We used data analysis parameters from the first inner- 
Galaxy observations of INTEGRAL'S core programme to esti- 
mate the possible precision of gamma-ray line centroid mea- 
surements. By using Gaussian fits to simulated spectra with 
varying levels of statistical noise, we tested the relation of the 
achievable line position determination accuracy and the signal- 
to-noise ratio. We confirm the centroid uncertainty to be pro- 
portional to the noise-to-signal ratio, therefore scaling with the 
inverse square root of exposure time. At present, noise is dom- 
inated by systematic uncertainties, resulting in a position un- 
certainty of +0.19 keV. If we extrapolate to the full exposure 
of 3 Ms scheduled for the first year of INTEGRAL opera- 
tions, we estimate a statistical uncertainty for the 26 A1 line cen- 
troid of +0.06 keV. We expect that systematic uncertainties in 
INTEGRAL results can be reduced in the near future through 
studies of SPI's spectral resolution, energy calibration stability, 
and background in-flight behaviour. When systematic uncer- 
tainties are reduced below the level of statistical uncertainty, 



we expect a 3<x detection of the Doppler shift relative to an 
equally long observation at the Galactic longitude of maximum 
Doppler shift at / = 330°. 

4. Conclusions 

We model 26 A1 sources in the Galaxy, adopting massive stars as 
the dominating sources. Using Galactic rotation and plausible 
assumptions for the fate of 26 A1 from ejection by the supernova 
or Wolf Rayet star into until decay in interstellar space, we de- 
rive an expected profile for the 1809 keV gamma-ray line from 
the decay of Galactic 26 A1. For this, we integrate along the lines 
of sight over the Doppler-shifted source regions of 26 A1 emis- 
sion with their intrinsic state of dynamical evolution. 

Measurements of the 1809 keV emission from galactic 26 A1 
play in principle a role similar to the Hi 21 cm radiation, with 
the benefit that the Galaxy is always optically thin to MeV 
gamma radiation. Current instruments are unable to improve 
the knowledge of galactic rotation over H i results, the basis fo r 
the inner galaxy rotation curve bvlOlling & Merrifieldif2000). 
but the knowledge of rotation can be used to test the position 
of 26 A1 sources. 

Our result demonstrates that the gamma-ray line profile re- 
flects the kinematics of decaying 26 A1 in the ISM. However, 
current Ge detectors will probably be unable to detect line 
shape departures from a simple Gaussian shape for spectra such 
as predicted by our model. Nevertheless, centroid and width of 
the line will provide a diagnostic for the 26 A1 sources in the 
Galaxy. 

Detection of the expected amount of line shift would put 
a limit on the contribution of local emission to the 26 A1 flux 
in the direction of the inner galaxy. By their massive star ori- 
gin, 26 A1 and H n-regions are connecte d; the latter having als o 
been used to measure galactic rotation (Brand & Blitz , 1993). 
H n regions are created by the ionising radiation from O stars, 
therefore probe the first 2 My of star formation. The peak of 
26 A1 emission from a coeval group of m assive stars only begin s 
at 2 My, continuing for another 10 My (Plusch ke et al.Ll200lT) . 
This demonstrates that 26 A1 and the other analyses of galactic 
rotation may complement each other. 

Our model parameters, chosen to represent supernova- 
produced 26 A1 with ejection into an typical ISM environment, 
predict a line width of » 1 keV for the inner Galaxy. This value 
is on the low side of the recent RHESSI measurement dSmithl 
2003), and significantly below th e "broad" line repo rted from 
the GRIS balloon measurement llNava et all I1996I) . We sug- 
gest, therefore, that 26 A1 decay in the interstellar medium is 
likely to teach us more about the ejection and expansion char- 
acteristics of nucleosynthesis ejecta, and thus about the mor- 
phology of the interstellar medium in the vicinity of massive 
stars, our assumed sources of 26 A1 in the Galaxy. For example, 
if 26 A1 sources typically are surrounded by cavities from earlier 
massive-star action, this could lead to some additional broaden- 
ing of the observed 26 A1 gamma-ray line. Further imaging spec- 
troscopy analysis of RHESSI data and of measurements with 
the S PI Ge spectrometer o n INTEGRAL (launched in October 
2002: IWinkler et alll2003l) promise to provide the data for such 
studies. 
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